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Preface 


The  investigation  reported  herein  was  conducted  for  the  US  Army  Engineer 
Waterways  Experiment  Station  (WES)  under  Contract  No.  DACW39-89-K- 
0018  by  Dr.  Alan  L.  Prasuhn,  South  Dakota  State  University,  Brookings,  SD, 
under  Work  Unit  32552,  "Sediment  Transport  in  Small  Channels,"  of  the  Flood 
Control  Channels  Research  Program  of  the  US  Army  Corps  of  Engineers  Civil 
Works  Research  and  Development  Program.  It  documents  modifications  to  the 
original  Ackers- White  sediment  transport  function  to  allow  for  multiple  grain 
sizes,  to  include  the  new  routine  in  HEC-6,  and  to  increase  the  number  of 
grain  size  dasses  into  the  cobble-boulder  range. 

The  study,  conducted  during  the  period  1989  to  1990,  was  under  the  gen¬ 
eral  supervision  of  Messrs.  F.  A.  Herrmann,  Jr.,  Director  of  the  Hydraulics 
Laboratory;  Mr.  R.  A.  Sager,  Assistant  Director  of  the  Hydraulics  Laboratory; 
Mr.  M.  B.  Boyd,  Chief  of  the  Waterways  Division,  Hydraulics  Laboratory;  and 
under  the  direct  supervision  of  Mr.  W.  A.  Thomas,  Research  Hydraulic  Engi¬ 
neer,  Waterways  Division.  This  report  was  prepared  by  Dr.  Prasuhn  as  part  of 
the  contract,  and  was  reviewed  by  Mr.  Thomas,  who  was  the  Contracting 
Officer’s  Representative. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


Conversion  Factors, 

Non-SI  to  SI  Units  of  Measure 


Non-SI  anils  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


1  Introduction  and 
Procedures 


The  puipose  of  this  report  is  to  indicate  what  has  been  accomplished 
concerning  die  modification  of  the  Ackers- White  sediment  transport  procedure 
and  its  incorporation  into  HEC-6.  The  original  proposed  changes  to  die 
Ackers-White  procedure  will  be  considered.  These  will  be  followed  be  a  dis¬ 
cussion  of  the  initial  problems  that  arose  when  die  modified  Ackers-White 
procedure  was  incorporated  into  HEC-6.  Finally,  the  procedure  that  evolved  as 
a  result  of  the  above  difficulties  and  recommendations  concerning  their  use  in 
HEC-6  will  be  presented. 

The  original  Ackers-White  procedure  (1973)  has  been  thoroughly  discussed 
in  die  literature,  e.g.,  Prasuhn,  Lewandowski,  and  Bagherzadeh  (1987).  The 
development,  other  than  that  necessary  to  explain  the  modifications,  will  not  be 
repeated  here. 

In  addition,  the  HEC-6  code  was  expanded  to  include  five  additional  sedi¬ 
ment  sizes  (from  small  cobbles  up  to  and  mrindiwg  large  boulders).  This  will 
be  covered  first  and  in  a  fairly  brief  fashion. 


2  Expansion  of  HEC-6  to 
Include  Five  Additional 
Sediment  Sizes 


The  existing  HEC-6  code  will  handle  15  size  fractions  of  sediment 
(U.S.  Army  Engineer  Hydrologic  Engineering  Center  1977).  The  research 
proposal  included  the  expansion  of  that  code  to  include  small  and  large 
cobbles.  Consistent  with  the  format  which  addressed  the  sand  and  gravel  sizes 
in  groupings  of  five  sizes,  the  expansion  actually  includes  three  additional 
sizes:  small  boulders,  medium  boulders,  and  large  boulders.  Subsequent  to 
that  report  a  complete  listing  of  all  the  required  changes  and  where  they  occur 
in  HEC-6,  subroutine  by  subroutine,  has  been  furnished.  They  were  also  high¬ 
lighted  in  a  copy  of  die  computer  program. 

The  expanded  program  now  accommodates  the  five  additional  sizes  and  die 
printout  should  be  consistent  with  all  previous  procedures.  The  dredging 
routines,  SDREDCE,  SSED,  JVOL,  and  perhaps  some  of  the  other  special 
options,  however,  were  not  thoroughly  tested.  There  should  be  no  operational 
problems.  The  print  format  should  be  checked,  however.  The  program  will 
now  handle  20  size  fractions  of  sediment  one  size  of  day,  four  sizes  of  silt, 
five  sizes  sand,  five  sizes  of  gravel,  small  and  large  edibles,  and  small, 
medium  and  large  boulders.  It  still  remains,  of  course,  to  find  a  transport 
function  for  the  large  sediments.  It  is  fdt  that  transport  of  at  least  die  cobble 
sizes  are  acceptable  using  this  procedure,  the  Schoklitscfa  procedure,  and 
maybe  the  Meyer-Peter  and  Muller  procedure. 

The  program  has  been  tested  with  many  combinations  of  sizes  up  to  the 
maximum.  This  induded  tests  with  both  four  and  less  than  four  silt  sizes. 
Changes  to  the  format  statements  have  been  avoided  with  few  exceptions.  The 
E-levd  printout  which  gives  all  the  Ackers-White  parameters  did  not  provide 
enough  room  for  large  values  of  D^.  This  column  has  been  changed,  but  die 
headings  (VF,  F,  etc.)  need  adjustment  too  as  they  are  inconsistent  with  the 
cobble  and  boulder  designations. 

Whereas  there  are  several  sediment  transport  functions  that  may  be  able  to 
give  a  reasonable  estimate  of  the  cobble  transport,  there  is  almost  no  verifica¬ 
tion  of  boulder  transport  If  HEC-6  is  used  for  the  transport  of  boulders,  care 

Chsplar2  Bqpsnaion  e4  HEC-6  to  Inetuda  Fh»  MMontl  Sedfcmnt 


should  be  taken  to  consider  the  reasonableness  of  die  results.  This  expansion 
of  the  code  wOl  not  be  discussed  further  in  this  report 
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3  Original  Ackers-White 
Procedure  (Summary) 


According  to  Ackers  and  White  (1973),  the  sediment  transport  depends 
upon  a  mobility  factor  Fp  given  by  the  following  expression 


Jgd^s-l) 


V 

\og(ay/ds) 


In 


(1) 


Here  s  is  the  specific  gravity  of  sediment  particles,  If*  =  the  shear  velocity,  g 
is  the  acceleration  of  gravity,  V  is  the  average  velocity,  y  is  the  depth,  ds  is  the 
representative  grain  diameter  (assumed  by  Ackers  and  White  to  be  the  *35 
sizes),  a  is  the  rough  turbulent  flow  coefficient  (assumed  by  Ackers  and  White 
to  equal  10),  and  n  is  a  factor  reflecting  the  sediment  size.  For  Dp  >  60, 
n  =  0,  otherwise  . 


n  *  1.0  -  0.56  log  Dp 


(2) 


The  dimensionless  grain  size.  Dp  is  given  by 


dp  -  <y«<»-iyv2]1/3 


p) 


where  v  is  the  kinematic  viscosity. 

Hie  dimensionless  transport  Gp  is  then  calculated  from  the  mobility  factor 
according  to 


Chapters  Original  Adwn-WMte  Procedure  (Summary) 


(4) 


Op  -  C(FplA  -  if 


where  the  coefficient  C,  the  exponent  m,  and  the  initiation  of  motion  parameter 
A  are  all  determined  by  regression  analysis  as  follows: 

Transition  range  (1  <  Dp  <  60) 


log  C  *  2.86  log  Dp  -  (log  Dp)2  -  3.53 

(5) 

m  -  9.66/Dp  *  1.34 

(6) 

A  =  0.23 /^Dp  *  0.14 

(7) 

Coarse  range  (Dp  >  60) 

C  *  0.025 

(8) 

m  *  1.50 

(9) 

A  -  0.17 

(10) 

The  resulting  dimensionless  transport  Gp  is  related  to  the  concentration  X 
(in  lb  sediment/lb  water,  or  N  sediment IN  water),  by 

Gp  «  (X/s)(y/ds)(u  JVf  (11) 


The  actual  sediment  transport  may  then  be  determined  by  either 

Gp^(A32)QfX  (tons/day)  (12) 

or 

Chapter  3  Original  Achere-WhUe  Procedure  (Summary) 


Ggr  -  &X 


(m 


(13) 


The  resulting  sediment  transport  is  consequently  the  total  bed  material  load 
based  on  the  representative  size.  It  was  not  suggested  by  Ackers  and  White 
that  their  procedure  could  be  used  to  calculate  the  transport  by  size  fractions. 


6 


Chapters  OrfQlrtel  Acfaara-Wilte  Prooadura  (Summary) 


4  Modifications  to  the 
Ackers-White  Procedure 
(Prasuhn  Procedure) 


The  original  proposal  to  modify  die  Adeem- White  proposal  has  been  pre¬ 
sented  and  discussed  on  several  occasions  (Prasuhn,  Lewandowski,  and 
Bagherzadeh  1987;  Heng  1989;  Prasuhn  and  Heng  1990).  In  its  final  form  it 
may  be  summarized  as  follows. 

In  their  procedure.  Adeem  and  White  introduced  an  initiation  of  motion 
parameter  A  which  they  expressed  as  a  function  of  the  dimensionless  grain 
size,  D-.  Using  primarily  transport  data  collected  by  White  and  Day  (1962) 
at  the  Wallingford  Research  Station  (designated  as  HRS  Series  A)  in  a  wide 
flume  with  a  broadly  graded  sample  a  similar  initiation  of  motion  parameter  A' 
has  been  determined.  The  A'  values  are  the  best  A  values  based  on  the  actual 
White  and  Day  measured  sediment  concentrations  followed  by  back  calcula¬ 
tions  for  A  using  die.  Ackers- Write  equations.  These  curves  are  plotted  in 
Figure  1  along  with  die  White  and  Day  data.  The  bed  distribution  of  these 
data  can  be  expressed  by  the  gradation  coefficient  a  =  (dg^/d^^2  «  4.198. 
The  original  Ackers  and  White  data  are  considered  to  work  best  for  a  single 
sediment  size;  thus  the  A  versus  D_  curve  is  associated  with  a  value  a  *  1. 
For  values  of  D„  <  56,  A'  is  greater  than  A,  reflecting  a  "hiding  factor" 
similar  to  that  of  Einstein  and  thereby  reducing  the  transport  of  the  finer  sizes. 
When  Dp  >  56,  A'  is  less  than  A,  resulting  in  an  "exposure  factor"  for  the 
larger  sizes  and  resulting  in  increased  transport. 

The  proposed  procedure  uses  the  a  value  of  a  given  bed  distribution  as  an 
interpolating  factor  to  get  a  best  A"  from  Figure  1,  for  each  size  of  bed 
material,  ami  hence  of  D_  For  each  grain  size  the  best  A"  is  calculated  as 
follows:  If  1  <  a  <  4. IS®,  then  A"  is  inteqxriated  from  die  respective  values 
of  A  and  A',  based  on  the  given  Z)_  If  a  «  1,  then  A'  «  A,  if  a  ■  4.196,  then 
A"  *  A'.  If  o  >  4.196,  then  A"  is  extrapolated  based  on  the  respective  values 
of  A  and  A'.  Thus  an  A*  is  determined  for  each  sediment  size.  Once  this  is 
accomplished,  the  transport  of  each  size  is  calculated  according  to  die 
conventional  Ackers- White  procedure.  The  individual  transport  is  then 
adjusted  by  the  pj  factor  to  reflect  the  actual  per  cent  of  that  size  found  in  die 
bed. 


(Praauhn  ftoo^uw) 


A  second  procedure,  to  be  referred  to  as  the  Heng  procedure,  involved  the 
use  of  two  graduation  coefficients,  Oj  and  defined  by  Oj  -  dj o/d16  and 
°2  =  ^ 94/d  jo-  These  were  used  as  interpolation  factors  in  the  same  way  as  a 
above.  They  were  used  to  reflect  the  bimodal  nature  of  many  bed  distributions 
and  in  tests  of  river  data,  did  a  slightly  better  job  than  the  Prasuhn  Procedure. 
However,  this  procedure  raised  additional  problems  when  incorporated  into 
HEC-6  and  further,  was  difficult  to  justify  theoretically.  Therefore,  although 
the  Heng  procedure  will  be  included  in  many  of  the  figures,  it  will  not  be 
recommended  nor  discussed  further. 
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5  Verification  of  the  Prasuhn 
Procedure 


Verification  of  the  modified  Ackers- White  procedure  has  beat  reported 
previously  (Prasuhn,  Lrwandowski,  and  Bagherzadeh  1987;  Heng  1989;  and 
Prasuhn  and  Heng  1990).  Additional  results  have  been  included  here  for 
selected  flume  and  river  data.  The  concentrations  shown  in  Figures  2-45 
represent  either  the  concentration  by  size  fractions  or  the  concentration  of  die 
total  bed  material  load  as  indicated.  There  is  clearly  a  variation  in  the 
accuracy  of  the  results,  but  they  are  generally  quite  acceptable  and  are  as  good 
as  can  be  achieved  by  other  transport  functions.  As  a  general  comment  the 
overestimation  of  the  very  fine  sand  and  the  tendency  for  gravel  sizes  to 
increase  as  the  particle  size  increases,  major  points  that  will  be  considered  in 
detail  lata,  are  not  obvious  here. 

Figures  2-14  represent  selected  Wallingford  HRS  Series  A  data  (White  and 
Day,  1982).  Since  the  A'  carve  was  based  on  these  data,  it  does  not  represent 
a  true  verification.  Figure  2  is  the  total  concentration  for  each  run  and  the 
remaining  figures  are  by  size  fractions.  On  these  figures,  and  those  that  fol¬ 
low,  the  L-B  results  can  be  ignored  as  representing  early  work.  The  A&W 
data  refer  to  the  original  Ackers-White  procedure,  and  the  Proposed  refers  to 
die  Heng  procedure.  Sediment  sizes  ranged  from  0.063,  die  beginning  of  size 
fraction  1,  up  to  15.7  mm,  die  end  of  size  fraction  8.  Those  results  labeled 
either  Prasuhn  or  A&W  will  usually  be  of  the  greatest  interest  Some  of  the 
Iowa  numbered  runs  refa  to  very  low  discharges,  but  generally  a  selection  of 
both  good  and  poor  results  will  be  included.  Note  that  each  pair  of  figures 
refers  to  a  specific  run,  in  one  case  referencing  the  calculated  to  the  measured 
data  by  a  line  of  perfect  agreement,  and  in  the  other  comparing  the  calculated 
to  the  measured  data  on  a  size  by  size  basis.  The  dashed  lines  in  the  first  case 
give  the  range  of  variation  within  a  factor  of  two.  The  size  by  size  compari¬ 
son  is  a  much  tougha  comparison  since  it  is  on  an  arithmetic  plot 

Figures  15-20  represent  selected  HRS  Series  B  (White  and  Day,  1982) 
results  for  which  the  sediment  sizes  ranged  from  0.063  to  635  mm.  These 
data  were  not  utilized  in  the  developmental  process.  Figure  15  is  once  again 
die  total  concentration,  but  based  on  the  calculation  by  size  fractions,  except 
for  the  original  Ackers- White  procedure  (A&W).  Note  that  the  original 
Ackers-White  procedure  overestimated  the  sediment  transport  in  both  cases. 
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Hie  remainder  of  the  figures  are  for  transport  by  size  fractions,  however,  ail  of 
the  runs  are  combined  in  Figure  16  which  pertains  to  just  the  Prasuhn  proce¬ 
dure.  In  some  cases  there  is  considerable  deviation  from  the  measured  data, 
but  not  in  any  consistent  fashion. 

The  remaining  flume  data  which  will  be  considered  here  were  the  data  sets 
collected  at  the  St  Anthony  Falls  Laboratory  and  reported  by  Hubbel  et  al. 
(1987).  The  generally  coarse  material  was  distributed  from  0.5  up  to  32  mm. 
These  results,  all  based  on  the  individual  size  fractions,  are  given  in  Fig¬ 
ures  21-32.  These  would  appear  to  give  an  excellent  verification  based  on 
independent  flume  data. 

Figures  33-42  are  based  on  sand  and  gravel  data  from  the  Platte  River 
(Kircher  1983).  With  the  exception  of  the  first  figure,  all  results  are  shown  by 
size  fractions  on  arithmetic  plots.  Both  die  Prasuhn  procedure  and  die  original 
Ackers- White  procedure  give  fair  results. 

Total  concentration  results  are  given  for  two  sites  on  the  Rio  Grande  River 
in  Figures  43  and  44.  Although  the  bed  material  was  mostly  medium  to 
coarse  sand,  the  sediment  ranged  from  very  fine  sand  to  medium  gravel 
(Nordin  1964). 

The  final  river  shown  (Figure  45)  is  the  Snake  River  data  (Jones  and  Sietz 
1980).  Since  the  measured  sediment  sizes  ranged  up  to  181  mm,  it  is  unfortu¬ 
nate  that  this  is  one  of  the  poorer  sets  of  total  concentration  results.  As  it 
turns  out,  the  original  Ackers- White  procedure  gives  the  best  results  which  is 
inconsistent  with  our  usual  observation  concerning  the  application  of  die 
Ackers- White  procedure  for  broadly  graded  sediments. 


^ - ■ — 
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6  Incorporation  of  the 
Prasuhn  Procedure  into 
HEC-6 


Problems  with  Incorporation 

The  original  process  of  incorporating  the  modified  Acker*- White  procedure 
into  HEC-6  involved  little  more  than  setting  op  a  separate  subroutine  which 
was  essentially  parallel  to  the  existing  Ackers- White  subroutine.  However  the 
incorporation  of  the  above  proposed  procedure  into  HEC-6  created  a  number 
of  problems  which  will  be  enumerated  here. 

Hie  problems  arose  during  either  the  coding  or  testing  stages.  They  will  be 
treated  separately  in  no  particular  order  below,  although  they  may  be 
interrelated.  The  resolution  of  the  problems,  or  recommendations  therefrom, 
will  be  considered  as  they  occur  as  well  as  summarized  at  the  end  of  die 
section. 


a.  The  gradation  coefficient,  a  *  (dg^d^^2  was  used  to  express  the 
spread  of  the  bed  sediment  distribution  and  as  an  interpolating  and 
extrapolating  parameter  for  the  incipient  motion  parameter  A*.  Previ¬ 
ously  only  sand  and  larger  sizes  woe  considered,  whereas  in  HEC-6 
significant  quantities  of  silt  and  clay  are  frequently  encountered.  It  was 
not  reasonable  to  include  the  clay  and  tilt  sizes  in  the  computation  of  a, 
so  Cor  purposes  of  this  computation  only,  the  bed  was  proportionally 
recomposed  tr  include  only  die  sand  and  larger  sizes.  This  can  not  be 
fully  justified  on  the  basis  of  hiding  and  exposure  factors,  but  no  alter¬ 
native  could  be  found.  It  does  not  appear  to  be  a  serious  problem. 

b.  In  the  computation  of  a  within  HEC-6,  evolving  bed  distributions, 
which  may  not  always  have  been  reasonable,  led  to  values  of  a  well 
beyond  what  had  been  previously  experienced  or  tested.  (The  HRS 
Series  A  data  has  a  value  o  *  4.198.)  It  was  felt  that  this  was 
nwTMiwrtir  and  o  was  set  equal  to  6  when  values  of  o  exceeded  6.  This 
again  seemed  to  work  wti*,  but  die  value  of  6  was  somewhat  arbitrary. 
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c.  During  the  testing  with  HEG-6  it  was  concluded  that  the  hiding  and 
exposure  effects  could  not  be  pronounced  when  only  sand  sizes  were 
present  so  the  modification  procedure  for  A *  was  bypassed  when  Z>  < 
SO  for  all  sizes.  This  increased  the  very  find  sand  transport,  making 
more  apparent  the  concerns  relative  to  the  overestimation  of  the  finer 
material  by  the  Ackers- White  procedure.  This  apparent  problem  will  be 
discussed  further,  below. 

d.  Some  considerable  evidence  is  now  available  concerning  the  overestima¬ 
tion  of  the  fine  sediments  by  the  original  Ackers- White  procedure.  In 
comparisons  based  on  the  HEC-6  test  runs,  both  the  original  Ackers- 
White  procedure  and  the  modified  Ackers- White  procedure  tended  to 
overestimate  the  transport  of  the  finer  material  and  very  fine  sand  in 
particular.  It  has  been  concluded  that  this  is  a  valid  criticism  of  the 
basic  Ackers- White  procedure  requiring  additional  analysis.  It  is  sug¬ 
gested  that  the  best  remedy  lies  in  the  adjustment  of  the  exponent  m. 

The  rationale  and  method  of  adjustment  is  Hkciimg**!  along  with  the 
coefficient  C,  below. 

e.  A  second  problem  coming  out  of  the  original  Ackers- White  procedure  is 
the  behavior  of  the  coarser  material,  primarily  gravel  and  larger  sizes. 
This  may  have  escaped  notice  previously  because  the  transport  was 
usually  limited  to  only  the  smaller  gravel  sizes.  However,  if  the  trans¬ 
port  is  normalized  so  that  each  size  is  considered  to  cover  100%  of  the 
bed,  transport  will  increase  as  the  sediment  size  increases. 

The  problem  which  is  presented  in  paragraph  e  is  demonstrated  in  Fig¬ 
ure  46  where  a  Froude  number  range  from  03  to  12  is  considered.  (The 
effect  becomes  much  more  pronounced  at  still  higher  Froude  numbers.)  A 
hypothetical  discharge,  and  width  were  picked  for  a  wide  rectangular  channel. 
TTre  required  depth  was  then  calculated  to  match  the  selected  Froude  numbers. 
Finally,  the  slope  was  determined  so  as  to  satisfy  the  Manning  equation  at 
constant  Manning  n.  The  calculations  for  sediment  transport  concentration  are 
based  on  the  assumption  that  each  sediment  size  (very  fine  sand  up  to  large 
cobbles)  completely  covers  the  bed. 

Except  for  the  smaller  Froude  numbers,  there  is  a  reversal  of  the  curve  for 
sediment  sizes  in  excess  of  Dp  equal  to  approximately  60.  Beyond  this  point, 
the  sediment  transport  increases  as  the  sediment  size  and  Dp  increase.  This, 
of  course,  is  entirely  unrealistic.  The  effect  is  more  pronounced  for  Fr  >  0.8 
(the  recommended  upper  limit  of  validity  as  recommended  by  Ackers  and 
White). 

Here,  there  is  little  to  go  on  except  the  illogical  behavior  of  the  current 
Ackers- White  procedure  with  regard  to  die  larger  sizes.  One  possible  fix  for 
this  problem  is  to  replace  die  existing  functions  for  C  with  alternative  equa¬ 
tions.  The  proposed  changes  to  C  will  also  be  discussed  below. 


sjnmptm  ©  lfw30Cp©f**Ofi  Of  WiO  nVMIWi  r 


Adjustments  to  m  and  C 


It  has  been  established  that  there  are  problems  with  both  the  Ackers- White 
procedure  and  the  modifications  thereto.  The  recommended  adjustment  for  m 
is  to  replace  the  existing  functions  with  an  alternative  expression  which 
reduces  the  magnitude  of  m  for  very  fine  sand  and,  to  a  lesser  degree,  fine 
sand  without  affecting  its  magnitude  for  the  larger  sediments.  The  following 
set  of  equations  were  tested  at  the  U.S.  Army  Engineer  Waterways  Experiment 
(WES)  to  determine  a  best  equation: 


m  -  8.419/D0,9  «■  1.289  (14a) 

m  *  8.027/D^0,9  ♦  1298  (14b) 

m  -  7.635/D  0-9  ♦  1.308  (l4*:) 

F 

m  -  7.244/D^0'9  ♦  1318  (14d) 

m  -  6.852/D^0,9  ♦  1328  (14e) 


These  equations  all  achieve  the  goal  of  reducing  the  transport  of  the  finer 
sediments  without  materially  affecting  the  transport  of  the  larger  sizes.  Equa¬ 
tion  14e  has  the  most  effect  In  each  case  the  exponent  of  0.9  was  chosen  to 
minimize  the  effect  on  the  larger  sizes. 

Typical  values  of  m  for  the  original  equation  (Equation  6)  as  well  as  Equa¬ 
tions  14a-e  are  given  for  various  values  of  Dp  in  Table  1. 

The  coefficient  C  was  likewise  to  be  tested  by  a  set  of  equations.  For  Dp 
values  up  to  166.7,  the  original  equation  (Equation  5)  was  used,  but  now 
extended  over  the  greater  range. 


t 

log  C  ■=  2.86  log  Dp  -  (log  Dp)2  -  333 


(5) 


Above  Dp  -  166.7  the  alternative  equations  were 
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0.00772 

(15a) 

03166/Dp  ♦  0.0QS82 

(15b) 

0.6296 /Dp  *  0.00394 

(15c) 

0.942 6/Dp  *  0.00207 

(15d) 

1235/Dp  *  0.000316 

(15c) 

Typical  values  of  C  for  D_  >  166.7  baaed  on  Equations  15a-e  are  given  for 
various  values  of  D_  in  Table  2.  Unfortunately,  adequate  time  was  not 
available  at  WES  tosatisfactorQy  test  all  of  the  above  equations  as  originally 
intended.  The  test  sets  that  were  completed  only  involved  the  sand  sizes  so 
the  range  of  equations  for  C  was  not  assessed.  The  only  tests  that  ran  success¬ 
fully  were  those  pertaining  to  the  HEG-6  Example  1  involving  the  Ozark 
Reservoir.  Example  1  ran  with  most  of  die  proposed  equations  for  m,  die 
original  Ackers-White  procedure  (MTC  «  7),  Test  Example  1  using  the 
Toffaled  procedure  (MTC  *  1),  and  Test  Example  3  which  is  identical  to  Test 
Example  1  except  that  die  Lanrsea  procedure  (MTC  *  3)  is  used  for  sediment 
transport.  Regrettably,  there  was  some  confusion  in  these  runs  snd  results  for 
only  four  of  die  five  m  equations  survive. 

The  IQ2C-6  computer  results  are  tabulated  in  ThMes  3-9.  Tribes  3, 4,  and 
5  are  the  Toffaled,  Laursen,  and  original  Ackers- White  remits,  respectively. 
Each  of  the  tables  contains  two  parts  reflecting  runs  of  1  and  64  days.  The 
original  Ackers- White  procedure  gives  higher  transport  of  dm  finer  material  as 
expected,  by  as  much  as  a  factor  of  2.  The  very  foie  sand  in  the  sediment  out¬ 
flow  perhaps  best  demonstrates  this  tendency.  The  effect  is  reduced  during  the 
second  time  interval,  but  this  involves  the  HEC-6  computations  and  iterations 
to  a  greater  extent  Tables  6, 7, 8,  and  9  refer  to  successive  equations  for  m. 
The  increased  reduction  in  very  fine  sand  is  apparent  in  each  table.  Again  the 
effect  is  reduced  for  the  second  time  period.  It  is  felt  that  dm  reduction 
afforded  by  Equation  14b  (Table  7)  gives  the  best  results,  and  this  becomes  dm 
recommended  equation  for  m,  although  additional  study  is  also  recommended. 

Although  the  tests  run  at  WES  did  not  provide  any  insight  into  the  behavior 
of  different  equations  for  C  on  HEC-6  output,  the  equations  had  been  tested  to 
some  extent  previously.  Equation  13e  'received  dm  greatest  attention,  and  on 
this  somewhat  limited  basis  Equation  15d  is  recommended  at  present.  This  is 
also  subject  to  the  recommendation  that  additional  testing  be  undertaken. 


The  above  changes  must  be  reasonable  and  consistent  with  the  original 
Ackers- White  procedure.  Refer  to  Figure  47,  which  is  a  copy  of  Ackers  and 
White’s  Figure  3  (1973).  Here  the  proposed  changes  to  m  and  C  have  been 
added  to  the  Ackers  and  White  graphs,  illustrating  that  the  proposed  changes 
are  consistent  with  their  original  data.  In  addition,  they  must  not  contradict  the 
soundness  or  logic  of  the  original  Ackers- White  procedure.  The  first  point 
was  somewhat  satisfied  by  testing.  The  second  can  be  justified  as  follows, 
based  on  the  development  by  Ackers  and  White  (1973).  The  numerical  values 
in  the  expressions  for  the  parameters  A,  n,  m  and  C  were  ultimately  determined 
by  regression  analysis.  They  (Ackers  and  White)  explain  that  the  regression 
first  led  to  the  expression  for  n  followed  by  the  incipient  motion  parameter  A. 
At  this  point,  the  expressions  for  m  and  C  were  determined.  Since  the  pro¬ 
posed  changes  in  m  and  C  are  in  mutually  exclusive  ranges  of  there 
should  be  no  violation  of  the  original  logic,  other  than  the  changes  in  m  and  C 
themselves. 
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7  Conclusions  and 
Recommendations 


No  nutter  which  version  or  form  of  Ackers- White  is  used,  the  argument  of 
die  log  function  in  the  expression  for  the  mobility  number  or  factor  (Equa¬ 
tion  1)  must  remain  greater  than  zero.  This  requires  that  y/ds  >  0.1.  This 
should  be  included  in  any  code  involving  Ackers- White.  For  example,  in  the 
present  subroutine  ACKER,  the  following  IF  statement  should  be  added  shortly 
after  statement  131  and  immediately  after  the  comment  line  "FOR  ■ 
SEDIMENT  MOBILITY  FACTOR11: 

IF  (EFD  LE  0  •  1*SD(I))THEN 
FGR*0 

GOTO  (die  write  statement  immediately  ahead  of  existing  statement 
number  statement  141) 

ENDDF 

Subsequently,  under  the  above  conditions,  die  calculation  for  GP(I)  needs  to  be 
bypassed  so  that  GP(I)  will  remain  at  die  initialized  value,  SPV.  (A  standard 
fix  up  in  the  existing  log  function  may  already  accomplish  the  same  goal,  but 
if  it  does  not,  this  should  avoid  problems.) 

The  expansion  of  the  code  to  cover  15  sand  and  larger  sediment  sizes  is  in 
operational  order.  It  should  be  useful  for  at  least  cobble  Led  streams  and  will 
increase  the  range  of  HEC-6  when  incorporated  into  die  program. 

Most  important  are  the  conclusions  and  recommendations  concerning  the 
Ackers-White  procedure.  It  is  still  felt  that  there  is  some  significant  gain  in 
the  modified  procedure  as  developed  at  South  Dakota  State  University.  How¬ 
ever,  it  is  now  felt  that  the  improvements  are  strictly  valid  only  when  the  bed 
distribution  is  broadly  graded.  It  may  still  be  desirable  to  incorporate  die 
program  into  HEC-6  under  these  conditions.  Although  reasonably  verified 
against  a  broad  range  of  river  data,  this  study  does  not  demonstrate  a  signif¬ 
icantly  improved  function  when  tested  in  HEC-6.  The  testing  did  identify 
problems  and  point  die  way  to  simple,  yet  direct,  improvements  to  the  Ackers- 
White  procedure.  By  replacing  die  functions  fin*  m  and  C  in  the  original 
Ackers-White  procedure,  much  of  die  criticism  of  the  procedure  is  eliminated. 
It  is  recommended  that  the  existing  subroutine  ACKER  be  used  in  HEC-6  with 


the  changes  made  to  the  equations  for  m  and  C  as  given  in  Equations  14b  and 
lSd.  If  the  modified  procedure  using  the  bed  coefficient  S  is  still  of  interest  to 
WES,  that  code  can  be  provided. 

By  the  aid  of  this  study  the  direction  of  work  had  changed  abruptly,  but 
for  the  better.  The  end  product  is  simpler  than  envisioned  at  the  start,  but  the 
recommended  changes  are  significant  Time  was  not  available  at  the  aid  to 
fully  evaluate  these  changes.  Consequently,  while  the  changes  will  improve 
the  performance  of  the  subroutine  ACKER  in  HEC-6,  they  may  not  yet  be 
optimized  It  is  recommended  that  the  equations  for  m  and  C  be  tested  further 


Chapter  7  Condudons  and 


References 


Ackers,  P.,  and  White,  W.  R.  1973.  *  Sediment  Transport;  New  Approach 
and  Analysis,*  Joum.  Hyd.  Div.,  ASCE ,  Vol  95,  No.  HY11,  pp  2041-2060. 

Heag,  H.  H.  1989.  'Sediment  Transport  Procedure  for  Nonuniform  Bed  Sedi¬ 
ment  by  Modification  of  the  Ackers  and  White  Procedure,*  Uresis,  South 
Dakota  State  University,  Brookings,  SD. 

Hubbd,  D.  W.,  et  al.  1987.  "Laboratory  Data  on  Coarse  Sediment  Transport 
for  Bed  Load  Sampler  Calibrations,”  USGS  Water  Supply  Paper  No.  2299, 
U.S.  Geological  Survey,  Denver,  CO. 

Jones,  M.  L.,  and  Sietz,  H.  R.  1960.  'Sediment  Transport  in  the  Snake  and 
Qearwater  Rivas  in  the  Vicinity  of  Lewiston,  Idaho,"  USGS  Open  File 
Report  80-690. 

Kircher,  J.  E.  1983.  "Interpretation  of  Sediment  Data  for  South  Plate  River 
in  Colorado  and  Nebraska,  and  Nodi  Platte  and  Platte  Rivers  in  Nebraska,” 
USGS  Professional  Papa  No.  1277-D. 

Nordin,  C  F.  1964.  "Aspects  of  Flow  Resistance  and  Sediment  Transport, 
Rio  Grande  Near  Bernalillo,  New  Mexico,”  USGS  Water  Supply  Papa 
No.  1498-H,  US  Geological  Survey,  Washington,  DC 

Prasuhn,  A  L.  1967.  Fundamentals  of  Hydraulic  Engineering,  Holt,  Rinehart 
and  Winston,  New  York. 

Prasuhn,  A.  L.,  and  Heag,  H.  H.  1990.  "Application  of  the  Ackers- White 
Procedure  for  Transport  by  Size  Fractions  and  Testing  in  HEC-6," 
Hydraulic  Engineering;  Proceedings  of  die  1990  National  Conference  on 
Hydraulic  Engineering,  Howard  H.  Chang  and  Joseph  C  HOI,  eds^  July 
30- August  3,  1990,  San  Diego,  CA,  Hydraulics  Division,  American  Society 
of  Civil  Engineers,  New  York,  pp  481-486. 

Prasuhn,  A  L.,  Lewandowski,  J.  A,  and  Bagherzadeh,  M.  H.  1967.  'Ackers- 
White  Equation:  Transport  by  Size  Fractions,”  Hydraulic  Engineering; 
Proceedings  of  die  1987  National  Conference  on  Hydraulic  Engineering, 


R.  M.  Ragan,  ed,  August  3-7,  1987,  Williamsburg,  VA,  Hydraulics 
Division,  American  Society  of  Civil  Engineers,  New  York,  pp.  936-941. 

U.S.  Army  Engineer  Hydrologic  Engineering  Center.  1977.  "HEC-6:  Scour 
and  Deposition  in  Rivers  and  Reservoirs,  Users  Manual,”  Davis,  CA. 

White,  W.  R.,  and  Day,  T.  J.  1982.  'Transport  of  Graded  Bed  Material,” 
Sediment  Transport  in  Gravel  Bed  Rivers,  R.  D.  Hey,  ed,  Wiley,  New 
York,  pp  181-223. 


Table  1 

Values  of  m  versus  Dar  for  Equation  14 


Eq.  number 


«  11.00 

6.17  2.31 

1.82 

1.58 

1.50 

14a 

9.71 

5.80 

2.35 

1.86 

1.59 

1.50 

14b 

9.33 

5.60 

2.31 

1.84 

1.56 

1.50 

14c 

8.94 

5.40 

zzr 

1.82 

1.58 

1.50 

|l44  8.56 

5.20 

2.23 

1.81 

1.58 

1.50 

14# 

8.18 

5.00 

2.19 

1.79 

1.58 

1.50 

Table  2 

Values  of  C  versus  Dnr  for  Equation  15 


1Ml7 


0.00772 


0.00772 


0.00772 


0.00772 


0.00772 


0.00772 


0.00645 


0.00279 


000 


0.00772 


0.00614 


0.00772 


0.00772 


0.00457 

0.00425 

0.00410 

0.00301 

0.00284 

0.00231 

0.00155 

0.00093 

0.00062 
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Figure  1.  Curve  for  A  vs 


Calc.  Concentration  (mg/l) 


Figure  3.  HRS  Series  A,  Run  S,  Concentration  by  Size  Fractions 


Figure  4.  HRS  Series  A,  Run  5,  Concentration  by  Size  Fractions 


Figure  6.  HRS  Series  A,  Run  7,  Concentration  by  Size  Fractions 
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Figure  7.  HRS  Series  A,  Run  8.  Concentration  by  Size  Fractions 


Figure  8.  HRS  Series  A,  Run  8,  Concentration  by  Size  Fractions 


Figure  10.  HRS  Series  A,  Run  9,  Concentration  by  Size  Fractions 


Figure  11.  HRS  Series  A,  Run  10,  Concentration  by  Size  Fractions 


Figure  12.  HRS  Series  A,  Run  10,  Concentration  by  Size  Fractions 


Figure  15.  HRS  Series  B,  Total  Concentration 


Figure  16.  HRS  Series  B,  Al  Runs,  Concentration  by  Size  Fractions 


1  2  3  4  5  6  7  8  9  10  1 1  12  13  14  15  16  17 

Size  Fractions 


Figure  17.  HRS  Series  B,  Run  3,  Concentration  by  Size  Fractions 


HRS  SERIES  B  RUN  5 
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Figure  18.  HRS  Series  B,  Run  5,  Concentration  by  Size  Fractions 
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Figure  19.  HRS  Series  B,  Run  6,  Concentration  by  Size  Fractions 


Figure  20.  HRS  Series  B,  Run  8,  Concentration  by  Size  Fractions 


Figure  22.  St  Anthony  Fate,  Run  1,  Concentration  by  Size  Freefone 


Figure  24.  St  Anthony  Fate,  Run  2,  Concentration  by  Size  Fractions 
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Figure  25.  St  Anthony  Falls,  Run  3,  Day  3,  Concentration  by  Size  Fractions 


Figure  26.  St  Anthony  Fafls,  Run  3,  Day  3,  Concentration  by  Size  Fractions 
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Figure  27.  St  Anthony  Falls,  Run  3,  Day  9,  Concentration  by  Size  Fractions 


Figure  28.  St  Anthony  FaRs,  Run  3,  Day  9,  Concentration  by  Size  Fractions 
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Figure  29.  St  Anthony  Fate,  Run  4,  Day  3,  Concentration  by  Size  Fractions 
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Figure  30.  St  Anthony  Falls,  Run  4,  Day  3,  Concentration  by  Size  Fractions 
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Figure  32.  St  Anthony  Fafls,  Run  4,  Day  6,  Concentration  by  Size  Fractions 


Figure  33.  Platte  River,  Total  Concentration 


Figure  34.  Platte  River,  Sample  1,  Concentration  by  Size  Fractions 
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Figure  36.  Rate  R href,  Sample  6,  CunuantwUon  by  Size  Fractions 
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Figure  37.  Platte  River,  Sample  8,  Concentration  by  Size  Fractions 


Figure  39.  Platte  River,  Sample  11,  Concentration  by  Size  Fractions 


Figure  40.  Platte  River,  Sample  12,  Concentration  by  Size  Fraction* 


Figure  41.  Ptatta  Rhrer,  Sampto  13.  Concentration  by  Sfcre  Fractions 
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Figure  42.  Ptatta  River,  Sample  14,  Concentration  by  Size  Fraction* 


Figure  43.  Rio  Grande  River,  Section  A-2,  Total  Concentration 


Figure  45.  Snake  River,  Total  Concentration 


Figure  46.  Concentration  Distribution  as  a  Function  of  the  Froude  Number 


Figure  47.  Ackers  and  White  (1973)  Figure  3 
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